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Abstract. We studied the relation between the distribution of sunspot groups and the Gleissberg cycle. As the magnetic field is 
related to the area of the sunspot groups, we used area-weighted sunspot group data. On the one hand, we confirm the previously 
reported long-term cyclic behaviour of the sum of the northern and southern sunspot group mean latitudes, although we found 
a somewhat longer period (P ~ 104 years). We introduced the difference between the ensemble average area of sunspot groups 
for the two hemispheres, which turns out to show similar behaviour. We also investigated a further aspect of the Gleissberg 
cycle where while in the 19th century the consecutive Schwabe cycles are sharply separated from each other, one century later 
the cycles overlap each other more and more. 
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1. Introduction 

The activity of the Sun has been studied in several ways for 
a long time. The oldest observed and recorded features on the 
solar disk are the sunspots. Regular sunspot observations were 
started by Galileo in 1610 and since then the data sets have been 
being continuously extended. Using these databases one has the 
opportu nity to study the stati stical properties of solar activity 
(see e.g. Viti nskv et al lll986"). The most prominent cycle in the 
sunspot series and in all solar activity is the 11 -year Schwabe 
cycle. The long-term trend in the amplitude of the Schwabe cy - 
cle is known as the secular Gleissberg cycle ( lGleissberd.1 1 945l) . 
which has a varying time scale of 80-120 years. 

This paper aims to study the relation between the distri- 
bution of sunspot groups and the Gleissberg cycle. Although 
the origin of the Gleissberg cycle is unclear, it is generally 
accepted that the interaction of the solar differential rotation 
and the magnetic field plays a basic role in the generation of 
all solar activity. From the observational data sets a correla- 
tion can be seen between secular variat ions in magnetic ac- 
tivity and differential rotation variations /Vitin skv et ail 1 1986t 
[Yoshimura & Kambry, 1993; Javaraiah, 2003). Although the 
role of the differential rotation in the long-term variation is 
not yet clear, several authors have investigated the interplay of 
differential rotation and the large-scale magnetic field, which 
drives the solar grand minima. Different modulational mech- 
anisms have been proposed in the recent literature, namely 
(1) fluctuations in the a-effect which could cause the dynamo 
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to stochastically flip and lead to modulation of magnetic en- 
ergy; (2) a dynamic or-effect forced modulation; (3) two-layer 
dynamos with meridional circulation; (4) nonlinear dynamics 
via macroscopic (e.g. Lorentz force) or microscopic (e.g. A- 
quenching) feedback (see a detailed review in Tobias, 20oJ). 

The following study may provide some additional observa- 
tional background for the investigation of solar dynamo models 
dealing with the description of the long-term evolution of solar 
activity. 

2. Data and method 

We need a detailed and homogenous sunspot catalogue to 
examine the fine structure of the long-term variation of the 
sunspot group distribution, so we use the sunspot data of the 
Greenwich Photoheliographic Results (henceforth GPR) such 
as the area, the heliographic position and the time of observa- 
tion. 

First, we would like to reanalyse the long term cyclic be- 
haviour of the (signed) sum of the mean latitudes of sunspot 
groups in the northern and southern hemispheres, respectively, 
but we weighted the latitudinal distribution of sunspot groups 
with their area. In a previous work, Pulkkinen et al. ( 1999) con- 
sidered 10-day averages of these sums and they found that 
the "magnetic equator" of the Sun displays a long-term vari- 
ation between northern and southern sunspot group latitudes 
with a period of P = 33900 + 950 days, and an amplitude 
of 1.31 + 0.13 degrees. They obtained similar results with the 
yearly averages. An important difference with respect to the 
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Fig. 1. The daily area- weighted latitude of sunspot groups as a 
function of time. 
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Fig. 2. The daily latitudinal area-weighted extension of activ- 
ity as a function of time. The vertical dashed lines denote the 
Schwabe minima. 



work of'Pulkkinen et al. 

1/1999!) 

is that, we use only Greenwich 
data. Consequently, our dat a ser ies is shorter (108 years) than 
of that of IPuUckinen et a (143 years), but is homoge- 

neous, so we can avoid the question of the compatibility of the 
measurements. 

It is well known that solar activity is correlated with the 
toroidal magnetix flux generated by the solar dynamo; basically 
the larger the area of sunspot groups, the greater the toroidal 
magnetic flux. Hence, in order to get a more realistic picture of 
the distribution of the toroidal magnetic flux, it is better to con- 
sider the latitudinal distribution of sunspot groups weighted by 
their area. For this reason, we define the following expression 
as the area-weighted mean latitude of sunspot groups: 



(A), = 



(1) 



where A, refers to the area of ;th sunspot group, O, denotes 
the latitude of ith sunspot group, A^, is the number of sunspot 
groups in the examined time intervals in day, and f refers to 
the length of the examined time interval. We consider 2-day, 
monthly and yearly (i.e. t - 2, t - 30, t - 365, respectively) 
area-weighted mean latitude data series for both the northern 
and southern hemispheres. On the one hand, we have not found 
significant diff'erences between the results obtained with differ- 
ent time intervals data sets. On the other hand, the average life- 
time of sunspot groups is about one month, so we concentrate 
only on the 30-day area-weighted averaged data. 

Secondly, we would like to investigate the correlation be- 
tween the Gleissberg cycle and the distribution of sunspot 
groups. In this case, we do not compose the time average, 



i.e. f = 1 . The daily area- weighted latitude of sunspot groups as 
a function of time is shown in Fig.^ This figure indicates that 
near the different Schwabe minima the distributions of sunspot 
groups deviate from each other. To analyse this discrepancy, 
we determine the daily latitudinal area-weighted extension of 
activity, which can be written in the foUowing form: 



i:fj,((A),-(D,)2 

Nr 



E, = 



(2) 



Fig. |3]illustrates the daily latitudinal area-weighted extension 
of activity as a function of time. In this case the discrepancy 
between the Schwabe minima is more significant, suggesting 
that as a Gleissberg cycle tends to its maximum, the consec- 
utive cycles overlap each other more and more, which will be 
discussed in the next section. 



3. Results and discussion 

3.1. Long-term behaviour of area-weighted mean 
latitudes of sunspot groups 

As it was mentioned in the previous section, we reanalysed the 
long-term variation of the (signed) sum of the mean latitudes of 
sunspot groups in the northern and southern hemispheres. An 
important difference to the work oflPuIkkinen et al. ( 199^ is 
that we considered the latitudinal distribution of sunspot groups 
weighted by their area. As a first step of our analysis we applied 
the discrete Fourier-transform (DFT) method for the data. The 
main peak belongs to the period of Pdft = 42 34 0+1 000 days. 
This is significantly larger than the period of Pulkkinen et alJ 
(1999), even larger than the length of our 39 000 day-long ob- 
serving window. As in such a case the fundamental frequency 
in the power spectrum is misplaced by aliases, we also searched 
for the period by phase-dispersion (PHD) method. This re- 
sulted in a period of Pphd - 37 900 + 100 days. The differ- 
ence between the two frequencies is |vdft - >'phdI = 3 x 10"^ 
day The reason can be understood from the properties of the 
spectral window of the data. The fundametal frequency of the 
DFT is vdft = 2.3618 x 10"^ day '. The full width at half 
magnitude of the peak, which depends on the length of the data 
series, is about 2 x 10"^ day"'. Consequently, the frequency 
found by the phase-dispersion method is located very close to 
the peak. Furthermore, as it is well known, if the spectral win- 
dow is symmetrical to the origin, a peak presents also at around 
-Vdft- This peak, as well as its side lobes can modify the exact 
location of the fundamental frequency in the spectrum. 

Supposing that the mathematical form of the sinusoidal 
variation is the following: 



/(f) - Go + cii cos(27rv0 + b\ sin(27rvf), 



(3) 



we list in Table[2the values of the Fourier coefficients obtained 
by least-squares fits, using the prev iously found periods, as well 
as the period given in the paper of lPuUckinen et alJ lll999l) . and 
the;^^ values obtained by test. 
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Fig. 3. The DFT-based power spectrum for the difference be- 
tween the ensemble average area of sunspot groups for the two 
hemispheres, i.e. (A3o/A^3o)North - (A3()/A^3o)south- 
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Fig. 4. Difference between the ensemble average area of 
sunspot groups for the two hemispheres, i.e. (A3o/A^3o)North - 
(A3o/A^3o)south- Each dot denotes a value for a 30-day inter- 
val, and the solid line represents the fitted sinusoidal pro- 
file prescribed in Eq. |3 with the following parameter val- 
ues: y = 1/48500, ao = 3.283199, ai = -6.971889 and 
bi = -18.475170. 

Table 1. The parameters of the sinusoidal curve fitted to the 
sum of the mean latitudes of sunspots of the two hemispheres 
by linear least squares for three different periods. Numbers in 
parenthesis denote the formal errors in the last digits. 



PHD DFT Pulkkinen et al. (1999) 

P[d] 37 900+ 100 42 340 ± 1000 33 900 ±950" 

ao n -0.331706(28) -0.383840(28) -0.303106(29) 

ai n 1.013833(39) 0.337501(40) 1414234(39) 

bi n -1.043614(39) -1.323759(39) -0.361780(40) 

y2 6.3164324 6.3197944 6.3222383 



" Note: For fitting we use only the period, because their fitted sinu- 
soidal profile did not contain a constant parameter (see Eq. 1 in 
IPulkkinen et alill999h 



3.2. Long-term behaviour of the ensemble average 
area of sunspot groups 

To confirm previous results we would like to investigate the 
long-term variation of the ensemble average area of sunspot 



groups. As we mentioned, the toroidal magnetic flux corre- 
lates to the sunspot area, so we construct another dataset in 
the following way. We take the total area (A30) and the num- 
ber of sunspot groups (A^3o) in every 30 days for the northern 
and southern hemispheres, respectively. We derived the differ- 
ence between the ensemble average area of sunspot groups 
for the two hemispheres, as (A3o/A^3o)North - (A3{)/A^3o)south- 
Applying the DFT and the phase dispersion method for these 
data, the same fundamental period was found (within the er- 
ror range), namely P - 48 500 + 100 days. The values of 
the Fourier-coefficients obtained by least-squares fits are aq - 
3.283199(29), fli = -6.971889(46) and /?! = -18.475170(36), 
while the x"^ value isx"^ = 130.0788751. 

Our results suggest that the toroidal magnetic flux shows a 
long-term variation between the northern and southern hemi- 
spheres, that is, if the ensemble average area of sunspot groups 
is larger on the southern hemisphere than on the northern one, 
then we can assume that the toroidal magnetic flux is also larger 
in the south, and later this dominance alternates. 

3.3. A further aspect of the Gleissberg cycle 

In Figs. lll2l it can be seen that the consecutive cycles overlap 
more and more around the Schwabe minima. In order to study 
the overlapping areas of the consecutive cycles on a long time 
scale, we investigated the behaviour of £, (see Eq. |5J around 
the minima of 11 -year cycles. Thus, we considered two-year 
intervals centered on the individual minima defined by t he rela - 
tive sunspot number (see review in Usoskin & Mursulal Eool . 
Then we calculated the mean and the standard deviation of the 
daily latitudinal area-weighted extension of activity ((£')min and 
C£,min, respectively) for these intervals. 

Fo r our study we invoke the results of lGarcia & MouradianI 
( Il998t) . They calculated the secularly smoothed sunspot num- 
ber of maxima R^fiN) and minima R*„ (N), where is the 1 1 -yr 
solar cycle number and they studied the variation of these com- 
puted values in time from the viewpoint of the Gleissberg cycle. 
Accordingly, their results are well comparable with ours. 

In Fig. |5j we present the mean ({E}jnin) and the standard 
deviation (cr^ „,;„) of the daily latitudinal area-weighted exten- 
sion of activi ty around the minirna of 1 1 -yr solar cycle and 
the results of lGarcia & MouradianI ( 1 1998|) . namely R*„(N) It is 
well visible that as a Gleissberg cycle tends to its maximum 
(£')min increases. While in the northern hemisphere this ten- 
dency is obvious, in the southern one it not so, but the tendency 
is visible. Both these numbers show that the consecutive cycles 
overlap each other more and more as the Gleissberg cycle ap- 
proaches its maximum, i.e. the sunspot groups of the next cycle 
appear before the end of the previous cycle. 

4. Conclusion 

In this paper we studied the long-term evolution of the distribu- 
tion of sunspot groups from different aspects, e.g. latitude, area 
and extension of activity. We applied two new methods, which 
in our opinion give a more realistic picture of the long-term 
variation of the toroidal magnetic field. Our method is based 
on the results of iPetrovav & Szakalv. (.1999.) . They found that 
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Fig. 5. The mean (solid line) and the standar d deviation (dashe d line) of the daily latitudinal extension of activity as a function of 
time. The dotted line represents the results of lGarcia & Mouradiaa(J998a) . i.e. the secularly smoothed sunspot number of minima 



the latitudinal distribution of the magnetic field at the surface 
reflects the conditions at the bottom of the convective zone, 
i.e. in this regard the convective zone behaves as a "steamy 
window". Although their original statement is restricted to the 
poloidal magnetic field, it can be assumed that this is also valid 
for the toroidal component (Petrovay, 2004, private communi- 
cation). This is why we applied an area-weighted method, since 
the sunspot area probably reflects the toroidal magnetic flux. 

First, we confirmed the results of IPulkkinen et al.1 ( 1 19991) 
using the area- weighted mean latitude data. However, we found 
a longer period applying the discrete Fourier-transform and 
phase-dispersion method, respectively. Based on the test, 
the period of Pphd - 37 900 days with an amplitude of 1 .04° 
turned out the best. We also studied the difference between the 
ensemble average area of sunspot groups for the two hemi- 
spheres, i.e. (A3o/A^3o)North - (A3o/A^3())south, which also indi- 
cates that the dominance of the toroidal magnetic flux oscil- 
lates between the southern and northern hemispheres in a long 
time scale. For this dataset we a found somewhat longer period 
{P - 48 500 days). Although the deviation of the two periods is 
greater than the formal errors (+100 days), the real uncertainty 
is significantly larger due to the shortness of the observing win- 
dow with respect to the length of the cycle and the large scatter 
Thus, we do not think that the difference of the periods implies 
real physical efi'ects. 

IPulkkinen et alJlll999l) interpreted this latitudinal variation 
as a mixed parity mode in which a quadrupolar component is 
oscillating with this period. They examined the toroidal com- 
ponents of both the quadrupolar an d the dipolar fields based 
on the paper of ^ Brooke et al.l lll998h . They claimed that mul- 
tiplying the absolute value of the total toroidal field by the 
signed mean latitude of the activity belt, this product should 
oscillate about zero. As the difference between the ensem- 
ble average area of sunspot groups for the two hemispheres, 
i.e. (A3()/A^3o)North - (A3o/A^3o)south is related to the strength and 
pattern of the toroidal magnetic field, our results may indi- 
cate observa tional evidence of such an oscillation. We note that 
IPipinI ( 1 19991) also found an asymmetry of magnetic activity. 

Beside studying the secular variation of the magnetic ac- 
tivity, we also investigated a further aspect of Gleissberg cy- 



cle. In Fig. lll2l it is visible that the sequential Schwabe cycles 
overlap each other at certain times. To characterize the rate of 
the overlapping we defined the daily latitudinal area-weighted 
extension of activity. We examined its mean and its standard 
deviation around the minima of the 1 1 -year solar cycle. These 
quantities are increasing as the Gleissberg cycle tends to its 
maximum (see in Fig. |5}. While in the northern hemisphere 
this tendency is obvious, in the southern one it is less so, but 
the tendency is visible. Nevertheless, as our investigation is 
limited to one Gleissberg cycle, we cannot declare unambigu- 
ous ly that this phenomenon correlates with the Gleissberg cy- 
cle. IPgl^^O {2000) found a similar tenden cy in the d i stance s 
in between nearby cycles (see Fig. 1 1 in IPelt et all 120001) . 
Hence, the continuous production of detailed, homogenous cat- 
alogues (e.g. Debrecen Photoheliographic Data successor of 
the GPR) is crucial from this point of view. On the other hand, 
the interpretation of the variable amount of overlap beetween 
the Schwabe cycle is a challenge for future solar dynamo mod- 
els. 
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